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Abstract 

The creation of large scale magnetic fields is studied in an inflationary universe where electrody- 
namics is assumed to be nonlinear. After inflation ends electrodynamics becomes linear and thus 
the description of reheating and the subsequent radiation dominated stage are unaltered. The 
nonlinear regime of electrodynamics is described by lagrangians having a power law dependence on 
one of the invariants of the electromagnetic field. It is found that there is a range of parameters 
for which primordial magnetic fields of cosmologically interesting strengths can be created. 



1 Introduction 

Magnetic fields are observed to be associated with most structures in the universe. Observations 
indicate magnetic fields on stellar upto supergalactic scales. The field strengths vary from a few //G 
on galactic scale, upto 10 3 G for solar type stars and upto 10 13 G for neutron stars. Furthermore, 
the magnetic field structure depends on the object it is associated with. Thus, e.g., magnetic fields 
observed in elliptical galaxies show a different structure from those associated with spiral galaxies 
&■ 

Magnetic fields in stars can be explained by the formation of protostars out of condensed 
interstellar matter which was pervaded by a pre-existing large scale magnetic field (see, e.g., [2]). 
An open problem remains to explain the origin of such large scale magnetic fields. 

There are different types of proposals. Ranging from processes on small scales, such as vortical 
perturbations and phase transitions to models taking advantage of the possibility of amplifying 
perturbations in the electromagnetic field during inflation in the early universe (see, e.g., [3]). 

Inflation provides a mechanism to amplify perturbations in some field to appreciable size. In 
order for this mechanism to lead to primordial magnetic seed fields of cosmologically interest- 
ing strength, the corresponding lagrangian should not be conformally invariant. The Maxwell 
lagrangian describing linear electrodynamics is conformally invariant. There have been already a 
multitude of proposals to break the conformal invariance of the Maxwell theory [3], e.g. by coupling 
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to a scalar field breaking Lorentz invariance [6], adding extra dimensions [7] or a coupling to 
curvature terms [8]. 

Here nonlinear electrodynamics is considered. It has its origins in the search for a classical 
singularity- free theory of the electron by Born and Infeld [9] . Later on it was realized that virtual 
electron pair creation induces a self-coupling of the electromagnetic field. For slowly varying, but 
arbitrarily strong electromagnetic fields the self-interaction energy was computed by Heisenberg 
and Euler (cf. PEj-|I2|). 

The propagation of a photon in an external electromagnetic field can be described effectively 
by the Heisenberg-Euler langrangian. Moreover, the transition amplitude for photon splitting in 
quantum electrodynamics is nonvanishing in this case. In principle, this might lead to observational 
effects, e.g., on the electromagnetic radiation coming from neutron stars which are known to have 
strong magnetic fields. \12\ I13j . In particular, certain features in the spectra of pulsars can be 
explained by photon splitting [13] . 

Finally, Born-Infeld type actions also appear as a low energy effective action of open strings 
[15\ [16] . As was shown in [T7] the low energy dynamics of D-branes is described by the Dirac-Born- 
Infeld action. 

The model of the cosmological background that will be considered consists of a stage of de Sitter 
inflation followed by reheating and a standard radiation dominated stage. Quantum fluctuations in 
the electromagnetic field are excited within the horizon during inflation. Once outside the horizon 
they become classical perturbations. As mentioned above, in general, the conformal invariance of 
the four dimensional Maxwell field has to be broken in order to amplify the perturbations in the 
electromagnetic field significantly. Thus, here electrodynamics is considered to be nonlinear during 
the de Sitter stage. This could be motivated by the presence of possible quantum corrections 
to quantum electrodynamics at high energies. However, once inflation ends electrodynamics is 
described by standard Maxwell electrodynamics. Thus the subsequent evolution described by the 
standard model of cosmology is unchanged. 



2 Nonlinear electrodynamics in the early universe 

The Born-Infeld or Heisenberg-Euler lagrangians provide particular examples of theories of non- 
linear electrodynamics. In general the action of nonlinear electrodynamics coupled minimally to 
gravity can be written as, see e.g. [HU HE] 

where L(X,Y) is the lagrangian of nonlinear electrodynamics. Furthermore, the invariants are 
denoted by X = \F^ U F^ V and Y = \F^ V *F^ u , where *F^ V is the dual bi- vector given by 
*F^ U = ^=e^ a PF al3 , and e^ ua/3 the Levi-Civita tensor with e i 2 3 = +1. 
The equations of motion are given by 

V M P^ = (2.2) 
where P^ u = — (LxF^ + Ly *F^ U ), furthermore La denotes La = dL/dA, and 

V M *F^ = 0, (2.3) 
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which implies that = d^A u — d v A il . The notation used is given by /i, v.. = 0...3 and = 
1, 2, 3. Moreover, the electromagnetic field is treated as a perturbation so that the vacuum Einstein 
equations apply to the background cosmology. The background metric is chosen to be of the form 



ds 2 = a 2 {r]) [-drj 2 + dx 2 ] . 



(2.4) 



Furthermore, following [4J the Maxwell tensor is written in terms of the electric and magnetic fields, 
E and B, respectively, as follows, 
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These are the components in the "lab" frame in which in linear electrodynamics the frozen- in 
magnetic field decays in an expanding universe as I /a 2 . 
Then equations (|2.2p and f|2.3|> imply 
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(2.6) 

(2.7) 
(2.8) 
(2.9) 



Although Maxwell's equations are recovered for lagrangians of the form F = n^X + mY, where no 
and m are constants, standard linear electrodynamics corresponds to F = —X. Equations (12. 6|) 
to (|2,9p is a set of four first order partial differential equations which can be transformed into a 
set of two second order partial differential equations. This describes the evolution of the nonlinear 
electromagnetic field in a curved background. In linear electrodynamics this procedure leads to 
two decoupled wave equations, one for the magnetic and one for the electric field. In the case of 
nonlinear electrodynamics the resulting wave equations for the electric and the magnetic fields are 
no longer decoupled because of the nonlinearities. 

Taking the curl of equation (|2.7p and using equations (|2.8|) and (|2.9[) a wave type equation for 
the magnetic field B can be found. 

±J?-fa*S) + 4^9,(a 2 B) + \^ dri {a 2 E) + ^ ( ^E - ^B 
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Similarly, taking the time derivative of of equation (|2.7|) and using the remaining equations results 
in a wave type equation for the electric field E, 



a 2 E) +8„ 
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A (a 2 E ) - d„ 
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(2.11) 



Equations (12. 1QH and (|2.1ip . respectively, are coupled and nonlinear which makes it quite diffi- 
cult to find exact solutions. However, as a first approximation it might be interesting to find the 
behaviour of the magnetic field neglecting the spatial dependence. This is the long wavelength 
approximation. Considering variations over a characteristic comoving length scale L much larger 
than the horizon aH then the spatial derivatives of a quantity can be neglected with respect to its 
time derivatives (see for example, [IH]). In general, E and B can be written in terms of Fourier 
expansions, 



E(x, rj) 



d 3 ke ik -* 



Ek(v) 



B(x, rj) 



B k ( V ). 



(2.12) 



Thus in the long wavelength approximation effectively only modes with small wave numbers will 
contribute to the Fourier expansions. Therefore, e.g. B(x,rj) ~ J^ c d 3 ke lk ' x Bk(rj). Just using one 



mode e %k ' x B^(r\) for k <C k c ~ aH one can show that in the limit k — > the terms involving spatial 
derivatives become sub leading. As it is commonly done, this approximation is applied to the second 
order equations (12. 1QH and (|2.11j) (cf., for example, [4]). 

A different way of looking at this is to use the stochastic approach to inflation where the mode 
expansion of a field, for example the inflaton, is separated into modes larger than the coarse graining 
domain and modes with wavelengths smaller than the coarse graining scale [20]. The superhorizon 
modes contribute to the coarse grained field which is made homogeneous by averaging over the 
coarse graining domain. The effect of modes leaving the coarse graining domain and contributing 
to the coarse grained field can effectively be modeled by a noise term in the equation of the coarse 
grained field. Neglecting this backreaction effect the dynamics of the field on superhorizon scales 
is basically described by the homogeneous, coarse grained field. 

Thus neglecting spatial derivatives equation (|2.10p implies, 



Lx Lx Lx 



Lx Lx 
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(2.13) 



where Bk 
r i i a t is ' 



a 2 Bk, £k = o?Ek and a prime denotes the derivative with respect to conformal time 
rj, tnat is ' = J^. In the case where the lagrangian only depends on X, Ly = 0, B k = const, is 
a solution which corresponds to the conformally invariant case, that is linear electrodynamics. In 
general, for Ly = 0, equation (|2.13p implies 



g, _Kk 
Bk ~ L X ' 



(2.14) 
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where Kk is a constant vector and Lx ^ 0. Moreover, linear electrodynamics is recovered for 
K k = 0. 

Furthermore, equation (|2,11|) implies 
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Equation (|2,15p can be integrated to give 

If _ £' _ 
L>X L>x 



(2.15) 



(2.16) 



where Pk is a constant vector. The homogeneous part of equations (|2.10p and (|2.1ip are coupled 
non trivially because of Ly, cf. equations (12. 13ft and (j2.15|) . Therefore in order to find solutions, 
the Lagrangian will be considered to be only a function of X, L = L(X). Furthermore, since 
X = \{B 2 — E 2 ) it is useful to find equations for £ k 2 and B% which are given by, for Pi > 0, 
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(2.17) 
(2.18) 



Assuming that the constant vector in equation (|2.16p vanishes, Pt = 0, leads to a significant 
simplification. In this case, equation (|2.16|) for L = L{X) can be solved immediately, giving for the 
electric field 



Mk 
Lx' 



(2.19) 



where Mk is a constant vector. Thus for Pk = equation (|2.18p leads to an equation only involving 
X and Lx, namely, 
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(2.20) 



In order to solve this equation a particular lagrangian has to be chosen. Here it is assumed that 
the lagrangian is of the form 



L 



x 2 \ — 



(2-21) 



where 5 is a dimensionless parameter and A a dimensional constant. This is the abelian Pagels- 
Tomboulis model |21j . The nonabelian theory was proposed as an effective model of low energy 
QCD [22]. Evidently, linear electrodynamics is recovered for the choice 8 = 1. The lagrangian 
(|2.2ip is chosen since it leads to a simplification of the equations, but still allows to study the 
effects of a strongly nonlinear theory of electrodynamics on the generation of primordial magnetic 
fields. In general, the energy-momentum tensor derived from a lagrangian L{X) is given by 



T =JL 

^ ~ 4tt L 



L x g a/3 F m F pv + 9ilv L 



(2.22) 
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Furthermore, for the lagrangian (|2.21|) the trace of the energy-momentum tensor is given by 



7T 



(2.23) 



which vanishes only in the case 5 = 1 that is for linear electrodynamics. In order to check if there 
are any constraints on the parameter 5 the energy-momentum tensor is calculated explicitly. The 
Maxwell tensor can be decomposed with respect to a fundamental observer with 4-velocity into 

an electric field E and a magnetic field B, following |23j, 



2E[ fM u u ] - r/^ K(7 u K B a , 



(2.24) 



where rj^ UKa = ^—ge^ Ufia and u^u^ = —1. Then the electric and magnetic field are given, re- 
spectively, by En = F l/tM u u and = ^v^k^^^"' The frame is defined by the proper lab 
coordinates (t, r) determined by dt = adij, df = adx. Applying a coordinate transformation then 
gives the relation between the fields measured by a fundamental observer and the lab frame. Thus 
using the four velocity of the fluid = (a" 1 , 0, 0, 0) this gives the relation [23] 



Ei — aEi, 



Bi = aBi 



(2.25) 



As shown in [23] the energy-momentum tensor of an electromagnetic field can be cast into the 
form of an imperfect fluid. The energy-momentum tensor of an imperfect fluid is given by (see for 
example, [23]), 



- /w 



(2.26) 



where p is the energy density, p the pressure, the heat flux vector and tt^u an anisotropic pressure 
contribution of the fluid. h^v = g^v + UpU v is the metric on the space-like hypersurfaces orthogonal 
to Up. With q^u^ = and ir^u^ = 0, 



r Ef J ,v'uf l h a 
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(2.27) 



Thus using equations (|2.22p and (12. 24ft the energy density and the heat flux vector for the Pagels- 
Tomboulis model (|2.2ip are found to be 
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Imposing the condition that ir a n is trace-free then the pressure and Tr a a are given by 
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(2.28) 
(2.29) 

(2.30) 
(2.31) 



Thus considering p (cf. equation (|2.28j) ) in general there is a constraint on 5. Namely, the positivity 
of p requires 5 > \. 
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Although in this work the abelian Pagels-Tomboulis model (cf. equation (|2.21|) ) is used, for 
completeness, other types of lagrangians describing theories of nonlinear electrodynamics are briefly 
summarized in the following. The self-interaction energy of a slowly varying, but arbitrarily strong 
electromagnetic field was calculated by Heisenberg and Euler |10l 111) . Expanding the resulting 
lagrangian into an asymptotic series gives a lagrangian of the form |10j-|12j 

L = X + k X 2 + k 1 Y 2 . (2.32) 

This describes the Heisenberg-Euler theory for the choice kq = and K\ = 45 Q 4 , where a 

is the fine structure constant and m e the electron mass. Assuming the coefficients to be general 
and imagining a situation where the quadratic term in X is dominant, the theory can be well 
approximated by the Pagels-Tomboulis lagrangian. 

Born-Infeld theory is another example of a theory of nonlinear electrodynamics. It was proposed 
as a classical, singularity-free theory of the electron [9]. The lagrangian is given by (cf. [15j-|17j) 

L = i (l - \A + W 2 X - /3 4 r 2 ) , (2.33) 

where (3 is a parameter. This type of action also appears in the description of open string states in 
string or M theory |15j-|17j. In this case [3 = 2-Tra' with a' the string tension. Considering a general 
parameter j3 and moreover the case in which the term [3 2 X is dominant results in a lagrangian of 
the Pagels-Tomboulis form. Furthermore, the parameter 5 in (|2.2ip is given by 5 = ^. 

However, since both invariants X and Y appear, the resulting equations are non trivially coupled 
in X and Y which makes it difficult to find solutions in closed form. In order to study the effects of 
nonlinear electrodynamics in the early universe, the abelian Pagels-Tomboulis theory (cf. equation 
(|2.2ip ) will be used. This has the advantage that even though it is a strongly nonlinear theory of 
electrodynamics, it is still possible to find approximate solutions in certain regimes. 

Therefore, in the following we will assume that the theory is determined by the abelian Pagels- 
Tomboulis lagrangian (|2.21[) . 



3 Estimating the magnetic field strength in the Pagels-Tomboulis 
model 

The following model will be considered. During de Sitter inflation electrodynamics is nonlinear 
and described by the Pagels-Tomboulis lagrangian (|2.21j) . This means that electrodynamics is 
highly nonlinear and very different from standard Maxwell electrodynamics. At the end of inflation 
electrodynamics becomes linear and thus the description of reheating and the subsequent radiation 
dominated stage is unaltered. 

The scale factor during de Sitter is given by 

a(??) = fll (£) \ (3-1) 

where rj < rji < 0. The end of inflation is assumed to be at rj = rji. 

Equations (|2.17|) and (|2.18|) are coupled, since X depends on E 2 and B 2 , in particular the invari- 
ant X reads, 2a A X ~ B 2 — ££. In order to make progress three different regimes of approximation 
will be considered. 
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i. ) ti 2 ~ o{£ 2 ). 

ii. ) B? ^> £u '■ This implies the approximation 2a 4 X ~ B?. 

iii. ) 3> B 2 . This implies the approximation 2a 4 X ~ — £ fc 2 . 

As will be explained further on, following [4 J it is assumed that quantum fluctuations in the 
electromagnetic field provide an initial magnetic and electric field. Therefore, it seems naturally 
to expect that initially ti\ ~ 0(£%). Thus the three cases mentioned above correspond to different 
types of evolution of the ratio B\j£\ during inflation. After the end of inflation, during the radiation 
era, while the electric field decays rapidly due to plasma effects, the magnetic field remains (see, 
e.g., HE5]). 



3.1 Case i.) ti>l ~ 0(8g) 

As a further simplification it is assumed that the contribution of the constant vector P£ in equation 
(I2.17P vanishes. In this case the equation (I2.20p during inflation in the Pagels-Tomboulis model 
(|2.2ip can be written in the form 
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where a dot denotes the derivative with repect to x and y 
definitions have been used, 



V 



Mi 



-l ' 



X 



y(x). Furthermore the following 

(3.3) 



K 2 



Moreover, the hatted quantities are dimensionless constants, 
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(3.4) 



Finally, Mp is the Planck mass. Equation (|3.2|) is a nonlinear differential equation and thus to 
find exact solutions is not trivial. Therefore, in order to proceed one further approximation will 
be made. It turns out that there is an approximate solution in closed form for 5 > 1. In this case 
neglecting the terms involving ( ^- j , with the exponents a = —4, —5, —6, yields the equation 



yy = $(y) 2 + 
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• my , 



where m 2 = — which is solved by 
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y(x) = C*2 [cosh (mx + (5 — l)mCx)} 1-5 , 



(3.5) 



(3.6) 
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where C\, C 2 are constants. In general, the magnetic field is given by 



M 2 / X \ - 2 ( S ~V 
% = £) ■ (3-7) 



5 2 a A VA 4 , 

However, for the approximate solution (|3.6p the first term becomes subdominant and the magnetic 
field can be approximated by 



x^j&Kjf) ■ (3 - 8) 

Thus, the magnetic field at the end of inflation at the time 771 can be expressed in terms of the 
magnetic field at the time, say 772, when the comoving length scale A was crossing the horizon 
during inflation. Thus 

B l( a i) e -w(\) cosh 2 [m(xi + (J- l)Ci)] 

B 2 (a 2 ) cosh 2 [m(rc 2 + ((5-l)Ci)]' 1 ' ' 

where N(X) is the number of e-folds before the end of inflation at which A left the horizon, that 
is, = a\jai. Furthermore, the constant C\ is chosen such that (5 — 1)C± = —x 2 - Using that 

during de Sitter inflation, a = 01(771/77) and the number of e-folds, results in the magnetic energy 
density ps at the end of inflation, 

Pb(ai) ^ PB(a 2 )e~ ANW cosh 2 [-mxi(e*W - 1)], (3.10) 

where pb = f^-- Following [1] the ratio of magnetic energy density to radiation energy density, r is 
introduced, 

r = ^-. (3.11) 

Pi 

In the case of linear electrodynamics, the energy density in the magnetic field and the radiation 
background decay both as a -4 and thus the ratio r stays invariant as the universe expands. In order 
to seed a galactic dynamo r has to be at least r = 10~ 37 corresponding to a magnetic seed field at 
the time of galaxy formation B s ~ 10~ 20 G. In order to seed the galactic field directly, without a 
galactic dynamo operating, r has to be of order r = 10~ 8 . Furthermore, we also note, that in a flat 
universe with a cosmological constant, these bounds can be lowered significantly. In this case, r has 
to be at least r = 1CP 57 to successfully seed a galactic dynamo [26]. This implies that the magnetic 
field at the time of galaxy formation has to be at least of order B s ~ 10~ 30 G. Following [4J, it is 
assumed that the energy density stored in the mode with comoving wavelength A is of the order 
of the energy density in a thermal bath at the Gibbons-Hawking temperature of de Sitter space. 

This leads to pb{o-2) — H 4 ~ (jM?) > where the constant energy density during inflation is given 

by M 4 . Finally, using that at the end of inflation the energy density in the radiation background 
is given by p 7 = M 2 T^ H , where Trh is the reheat temperature, the ratio of magnetic to radiation 
energy density at the end of inflation at 77 = 771 is found to be, 



r(a 1 ) 



' M ' 
,M Pj 



~ 2 e -m» cosh 2 \- mXl (e"W - l)] . (3.12) 
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Furthermore, following [3], the number of e- folds can be found as 



~ 9.2 x 10 25 



Mpc 



M 



\Mp) 



(3.13) 



where it is assumed that the scale factor today is ao = 1 and thus comoving and physical scales 
coincide in the present. Thus the ratio r at the end of inflation is given by 



r{a\) ~ 10 
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Therefore, in order to achieve, a ratio of magnetic to radiation energy density at the end of inflation, 
which is at least some value ro, that is r(a{) > rg, mx\ has to satisfy, 



mx\ > 10 



-26 



A 



Mpc 
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(3.15) 



Equation (13.15h only implies a constraint on mx\ if the argument of arccosh is bigger or equal to 
one which can be interpreted as a bound on the reheat temperature. Assuming a galactic scale, 
that is A = IMpc implies, 



10 



52 



TrhX 
M J 



> 1, 



(3.16) 



which for ro = 10~ 37 implies Tpjj > 1 MeV, where it was assumed [1] that the inflationary energy 
scale is given by M = 10 17 GeV. This is always satisfied since the reheat temperature has to be at 
least 10 MeV (for even lower values, see |27| ) in order to allow for nucleosynthesis to take place 
unaltered. For a smaller value of ro, say ro = 10 -57 , the reheat temperature is required to be at 
least of order 10 3 GeV. However, in this work the reheat temperature is assumed to be at least 
T RH > 10 9 GeV. There is also an upper bound on — mx\ coming from the requirement that r < 1 
in order not to overdose the universe. This implies, 



— mx\ < 10 
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\MpJ 



arccosh 



10 
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(3.17) 



The constraint equations (13, 15ft and A3. lTf> can always be satisfied, since, by assumption, ro < 1 
and, moreover, for physically interesting models ro <C 1. Thus assuming A = 1 Mpc and M = 10 17 
GeV, the following values for — mx\ are found. For a model with reheat temperature Trh = 10 9 
GeV [3] the parameter — mx\ has to be in the range 2.7 x 10~ 20 < —mx\ < 5 x 10 -20 in order to 
achieve a magnetic seed field with a field strength to be at least B s ~ 10 _20 G, corresponding to 
ro = 10 -37 . For a higher reheat temperature Trh = 10 17 GeV [[4], for the same magnetic seed field 
strength —mx\ has to be in the range, 9.5 x 10 -23 < —mx\ < 1.5 x 10 -22 . And similarly, for the 



less conservative bound ro = 10 , for Trh = 10 GeV, — mx\ has to be in the range 1.4 x 10 < 
-mxi < 5 x 10" 20 and for T RH = 10 17 GeV it is found that 6.7 x 10~ 23 < -mx\ < 1.5 x 10~ 22 . 

Thus there is a range of parameters for which strong enough magnetic seed fields can be created 
in the Pagels-Tomboulis model of nonlinear electrodynamics. Since the analysis is based on the 
approximate exact solution given by equation (|3.6j) it is also important to check that the solution 
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provides a good approximation to the solution of the full differential equation (|3,2p . This has been 
done in Appendix A. 

In summary, for 5 > 1, there is an approximate analytical solution which allows to find an 
expression for the ratio of the magnetic to radiation energy density. There is a range of parameters 
for which magnetic seed fields of cosmologically interesting field strengths can be created. 

3.2 Case ii.) B% > S 2 

In this case equation (12. 18ft takes the form, 

!^*> + <'-4£^-£(£) = * (3 ' 18) 

where it has been used that 2X ~ B\. It is possible to find different types of solutions of equation 
(|3.18p depending on the value of the parameter 5 of the Pagels-Tomboulis model. On the one hand 
there are power law solutions for 5 ^ \ an d S ^ |. On the other hand there are solutions with 
a distinct behaviour for 5 = ^ and S = |. Actually of the latter ones only the case 5 = \ will 
be discussed explicitly. This is so because for 5 = § it is only possible to find an implicit solution 
depending on Euler's (3 function which makes it very difficult to estimate the primordial magnetic 
field strength. 

3.2.1 Solution for 6 ^ \ and S ^ f 

For 5 7^ \ and 5 ^ |, equation (|3.18h can be solved by a power-law function, 

(3-19) 

This leads to the solution for the magnetic field 



B 2 k = 2A 4 



i 

5 - 4<5 N 2 



02^1 \2(5 — 1 



1-2,5 / x 157 — j- 



m 



(3.20) 



where, as before, a 2 = J^4 S 2 ■ Thus using the definitions as given for case i.) (cf. section |3~T|) the 

ratio of magnetic energy density to radiation energy density at the end of inflation r(ai) is found 
to be, for S ^ ±, 8 ^ f , 

The range of validity of the assumption B£ 3> £ 2 can be checked to first order by using the 
solution for B\ ~ 2X, (cf. equation (|3,20p ) in the equation for £ fc 2 , equation (|2.17p . This leads to 

4 2 " + lV -i£Z ' + f3 2V ~% 2 = 2Pl (3.22) 

where (3\ = 3a(5 — 1) and 02 = 2a{5 — 1) [a (6 — 1) — 1]. This equation is solved by, for P% ^ 0, 

9 -4.5+5 12(^-1) 



ii 



where Co and c\ are constants. For P| = 0, the solution is given by equation (|2.19p which will be 
discussed below. During de Sitter inflation the scale factor is given by equation (|3.ip . Thus, finally, 
the ratio £j?/B£ is given by 

where fio, [i\ and [12 are constants depending on the constants of the solutions of the electric 
and magnetic field. However, their explicit form is not important here. Since during inflation, 
77 < 771 < 0, and hence 77/771 > 1. At 77 = 772, that is at the time when the mode is leaving the 

horizon during inflation, the initial condition, -^-(772) = 1 is imposed. In order for the approximation 
to be consistent, it is required that the solution evolves such that < 1. So assuming that each 

"k 

of the terms is of the order of | at 77 = 772 then the constants fio, \i\ and 7x2 can be estimated in 
terms of — . This leads to 



% = l(n) U ** + l(lL) m + \(lL ] . (3 . 25) 



Thus the last term is growing and in general, is not smaller than one. Thus in order for 

"* 

the solution to be consistent within the approximation, the constant c\ in equation ([3.23P has to 
be set to zero. Furthermore, the exponents in equation (|3.25p have to be positive, imposing the 
constraints, 5 < \ or 5 > f . The former one is ruled out since 5 > \. Then by assuming that the 

two remaining terms, after setting c\ to zero, contribute equally at 77 = 772, the ratio < 1 for 

7/ > 772 and 8 > § . 

In order to seed the galactic dynamo, it is required that r > r$ where ro is the lower bound 
on the strength of the magnetic field. In the expression for r = ^ at the end of inflation there 
are, apart from 5, two parameters: on the one hand the constant energy density during inflation 
given by M 4 and on the other hand the reheat temperature, Tjoj- Following [1] M is chosen to 
be M = 10 17 GeV. The reheat temperature depends on the details of the reheating process. It 
can be as low as 4 MeV [27] and in general, one expects an upper limit T^h < M. However, in 
supersymmetric theories this limit is lowered down to 10 9 GeV [28j . 

In Figure 1 log r is plotted against the Pagels-Tomboulis parameter 5 for the inflationary energy 
scale M = 10 17 GeV for the reheat temperatures Trh = 10 17 GeV and Trh = 10 9 GeV. As can 
be seen from Figure 1 for S > | there is range of S for which primordial magnetic fields with 
cosmologically interesting field strengths can be generated. 

In the case P? = the solution for £^ is given by equation (|2.19p . This leads to 

'! * it iw) (3 - 26) 

and thus 

4 2 .. 1 (s-*sY{vY 2 (3 . 27) 



Bl m 2 x\ \25 -lj V 7 /! 



12 



-20 



-40 



-60 



-100 









/ / 




1 / 





1.5 2 2.5 3 3.5 4 4.5 5 

6 



Figure 1: For A = 1 Mpc logr (cf. equation (|3.2ip ) is shown as a function of 6 for Trh = 10 17 GeV 
(black line) and Trh = 10 9 GeV (long-dashed line). The dashed line corresponds to r = 10~ 37 . 

where as before (cf. equation (|3.2p ) m 2 = sa. Thus imposing the initial condition ^§-(772) — 1 
implies 



(3.28) 



which implies d| > 1 for 77 > 172- Thus this solution is not consistent with the approximation. 



3.2.2 Solution for 5 = ± 



In this case equation (|3.18j) can be written as 

1 ^ , g 



(3.29) 



where g = (jf\ X, x = rj/Mp 1 and 02 a constant as defined before, (cf. equation (|3.3|) ). This 



is solved by g(x) = C2 cosh 2 
B\ ~ 2X it follows that 



, where c\ and C2 are constants. Thus with 



Dz - 2c 2 ( ^ ) n.slr 



1 / \ 3 

is) X1 UJ +C1 



(3.30) 



This then leads to the ratio of magnetic energy density to radiation background energy density 



r = — at the end of inflation, 



-mi / A VV M V cogh 2 



r(ai) ~ 10- iU4 — — 



x 10 xi( — 1 



' A \ 3 / M \ 2 T, 



18/ V M P C / \Mp) M P 



RH 



. (3.31) 
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Here the constant c\ has been chosen as c\ 
ro < r(ai) < 1 results, for A = 1 Mpc, in 



10" 



M 
W P 



fT*g\ 
\M P J 



-i 



arccosh 



f^LY 2 ( Trh \ 



10 



52 



5 

3 I 
2 




< 10 



-78 



\Mp) \Mp) 



m 

-i 

arccosh 



Imposing the condition that 



< — Xi 



10 



52 



(£) 



(3.32) 



The resulting values for different choices of Trh, M and ro are given in Table 1. Furthermore, 



T RH (GeV) 


M(GeV) 




Xl ( ll)low 


Xl ( ii)u P 


10 9 


10 17 


io- 37 


8.6 x 10~ 63 


1.6 x 10" 62 


10 9 


10 17 


lO" 57 


4.4 x 10" 63 


1.6 x 10~ 62 


10 17 


10 17 


IO" 37 


8.6 x 10~ 71 


1.6 x 10" 70 


10 17 


10 17 


10~ 57 


4.4 x 10" 71 


1.6 x 10" 70 



Table 1: Lower and upper bounds of —x\ (-[f) 2 derived from equation (|3,32p for different values 
of the reheat temperature Trh, the constant energy density during inflation determined by M and 
the lower limit on the field strength of a primordial magnetic seed field determined by ro. The 

notation used indicates — x\ (ff) ; 2 < —x\ (y|) 2 < — x\ (y|) ; 



i 

' up 



in order to check the validity of the solution (|3,30p which was derived under the assumption that 
£ fc 2 /£> 2 <C 1, we consider the cases P| > and P% = 0. In the case P? > the electric field 
strength is determined by equation (|2,17p . As a first approximation, the solution for X ~ \B\, 
where -B 2 is given by (|3.30p . will be used in (|2.17p . For consistency, the resulting solution for the 
electric field strength should be much smaller than the magnetic field strength. In equation (|2,17p 

L' L" 

the expressions for and are needed which are given in Appendix B. The cosmologically 

interesting values of // = — (f§) 2 x\ are very small, /j < 0(1O~ 62 ), as can be seen from Table 1. 
Thus to zeroth order in fi equation (|2.17p becomes, 



c 2 // 



6 P 2 / , 12 ? 2 



(3.33) 



which is solved by 



P 2 n 2 



(3.34) 



where (5q and (5\ are constants. Therefore 

p2 2 / v_ 



4 2 



2c2d\ cosh 



Xl 



(3.35) 
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E 2 £ 2 
Imposing the initial condition -jjk\r)2) = 1 it can be seen that -J^ is decreasing very fast and thus 

"k "k 

the solution (|3.30p is consistent at lowest order in [i. 

In the case P% = the solution for the electric field strength is given by equation (|2.19p , leading 

to 



2X1 



(3.36) 



This implies 



5? 



T 



(3.37) 



Imposing that initially ^-(772) = 1 it follows that 



B 2 k 



(3.38) 



which is smaller than one for r/ > r/2- Hence the solution is consistent with the approximation. 

Thus, for the solution in this case, there is a choice of parameters which allows to create 
primordial magnetic fields of cosmologically interesting field strengths. This holds for both cases, 
P 2 > and P 2 = 0. 

3.3 Case Hi.) E 2 > B fc 2 

Starting out at the same order of magnitude, in this approximation at the end of inflation the 
energy density in the electric field is much larger than in the magnetic field. The approximation 
implies that £? ~ — 2a 4 X. In the case P 2 > equation (|2.17p yields to 



dr] 2 v ' ' X dr) 



(a 4 X) + 3(5 - 1)^4- (« 4 ^) + 2 ( 5 " 1) 



a A X = -Pi, (3.39) 



which is solved by, for 5^1, 



p2 2 



X 



2a\{Q5-b) 2 \rn 



(3.40) 



Thus using X = X\ [J}[j i n the equation determining the magnetic field (cf. equation ([2. 18ft ) 
gives, 



Bl 



a\5 2 [l-a(5-l)f \A 8 
where bo and b\ are constants and X\ 



V 2\ -(S-l) / \ 6-2a(<5-l) 

+ »0 



Vi. 

phi 



4 / \ 5-a(5-l) 

7] N 



24(65-5)' 



+ 6l 

ViJ \Vi 



-. With a = 6 this leads to 



, (3.41) 



-12(5-1) 



+ M1 — 



/'2 ( — 



5-65 



(3.42) 
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where fiQ, fj.% and \i2 are constants which can be found from the expressions for E? and -Bj?. Imposing 
the initial condition ££(772) — B?(j]2) and that all terms contribute equally at this time results in 



El 



12(5-1) 



3 V V 



3 \ r] 



6<5-5 



(3.43) 



Thus in order to achieve, ByE^ < 1 the constant bo in equation (I3.4ip has to be set to zero. With 
the remaining two terms contributing equally at 77 = 772 and requiring ^ < 5 < | leads to solutions 
which are consistent with the assumption B^/Ef, < 1. Furthermore, in the expression for B% the 
dominant contribution comes from the last term, thus the evolution of the magnetic field is given 



by B\ 



where (3 = 11 



65 and ^ < 5 < | . Moreover, the ratio of the energy density in the 



'k ~ \% 

magnetic field and the background radiation r at the end of inflation can be calculated, resulting 
in 



r(ai) ~ (9.2 x 10 25 ) 



Mpc 



- p ( M V 
\M~p) 



2j£ 
3 



TrhX 
Mp 



(3.44) 



In Figure 2 log r is shown. As can be seen the resulting magnetic field strengths are far below the 
lower boundary of ro = 10~ 37 , corresponding to a mag netic seed field of B s = 10" 20 G. 
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Figure 2: For A = 1 Mpc logr (cf. equation (|3.44p ) is shown as a function of 5 for Trh = 10 17 GeV 
(black line) and Trjj = 10 9 GeV (long-dashed line). The dashed line corresponds to r = 1CP 37 . 
The left panel corresponds to the case P 2 > 0. The right panel corresponds to the case P 2 = 0. 

Finally, the solution for P 2 = will be discussed. Thus using equation (I2TT9]) and X ~ -\E\ 
yields to 




(3.45) 



It is found that the solutions are consistent with the assumption > B% for 1 < S < | implying 
f3 = 4 and for 5 > | corresponding to (3 = 2||^j. Thus using (3 = 4 and A = 1 Mpc yields to 
r(ai) ~ 10" 104 for M = 10 17 GeV and T RH = 10 17 GeV. Moreover r(a x ) ~ 10" 77 is found for 
M = 10 17 GeV and T RH = 10 9 GeV. These values are far below the lower bounds on the primordial 
magnetic field required to seed the galactic field. The results for S > | are shown in Figure 2. As 
can be seen for Trh = 10 9 GeV magnetic fields satisfying r > 10~ 37 can be generated for 5 > 19.5. 
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3.4 Discussion 

Solutions for the magnetic energy density of nonlinear electrodynamics with a lagrangian given 



3-1 



by L = - 




X, where A and 5 are constant parameters have been found for different 



approximations. The solutions are determined by the system of equations (I2.17|) and (j2. 18|) . These 
equations depend on two constants, Pj: and a|. In the case where P% = 0, equation (|2.17|) 

_ 2 

is replaced by equation (12.191) which involves a new constant, , in the final equation (|2.20p . 
Furthermore, M| and K 2 lead to the definitions of the two dimensionless constants a% and oli in 
equation (|3.2p and m? = 

It is assumed that the electric and magnetic fields, respectively, have their origin in quantum 
fluctuations during inflation. Therefore, it seems quite natural to impose that initially, that is at the 
time when the perturbation leaves the horizon, the energy density in the electric and magnetic field 
are of the same order. During the later evolution these quantites of course can be very different. 
In order to solve the equations, we have assumed three different types of evolution of the ratio of 
the energy densities in the electric and magnetic field. This has led to estimates of the primordial 
magnetic field at the time of galaxy formation. 

In the case B? — 0{££) for P£ = it was found that strong primordial magnetic fields can be 
generated. 

Assuming that during inflation B 2 3> £ fc 2 there is a range of the Pagels-Tomboulis parameter 6 
for which in the case Pj 2 > primordial magnetic fields can be generated that are strong enough 
to seed the galactic dynamo. In particular, for 6 > 1.9 for Trh = 10 17 GeV and for 5 > 3.0 for 
Trh = 10 9 GeV the ratio of the energy density of the magnetic field over the energy density of the 
background radiation r is found to be r > 10™ 37 corresponding to a primordial magnetic field of at 
least B s = 10 _20 G (cf. Figure 1). However, in the case P 2 = this solution is not consistent with 
the approximation Bj 2 3> S 2 . Thus it cannot be used to estimate the primordial magnetic field in 
this case. The former class of solutions do not inlude the case 5 = ^. In that case the solutions 
found for the electric and magnetic field are consistent with the approximation for P? > and 
P 2 = 0. Moreover, the resulting magnetic field is strong enough to seed the galactic dynamo. 

Finally, making the approximation £ fc 2 3> B^ yields in the case Pj 2 > to very weak magnetic 
fields. However, in the case P| = 0, for 5 > 19.5 and a reheat temperature Trh = 10 9 GeV 
primordial magnetic fields result which could successfully act as seed fields for the galactic dynamo 
(cf. Figure 2). 

4 Conclusions 

Observations of magnetic fields on large scales provide an intriguing problem. A possible class 
of mechanisms to create such fields is provided by inflationary models. Fluctuations in the elec- 
tromagnetic field are amplified during inflation and provide a seed magnetic field at the time of 
structure formation which might be further amplified by a dynamo process. In general a sufficiently 
strong initial field strength can only be achieved if the conformal invariance of electrodynamics is 
broken. This has been realized, for example, in models where the Maxwell lagrangian has been 
coupled to a scalar field, to curvature terms, etc. or by breaking Lorentz invariance or adding extra 
dimensions. 

Here nonlinear electrodynamics has been considered. It has been assumed that whereas dur- 
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ing the early universe electrodynamics is nonlinear it becomes linear at the end of inflation. In 
particular the evolution of the magnetic energy density has been studied in a model of nonlinear 



s-i 

\2 



2 Tp rpfiu 
1 iiv L ) 



electrodynamics which is described by a lagrangian of the form L (F^F^) 2 /A 8 

where A and 5 are parameters. Originally the nonabelian version of this model had been proposed 
to describe low energy QCD [22j . Here this model has been chosen as it is a strongly nonlinear 
theory of electrodynamics which allows to study in a semi-analytical way the possible creation and 
amplification of a primordial magnetic field during de Sitter inflation. This is so since on the one 
hand the lagrangian only depends on one of the electromagnetic invariants, namely X = jF^F^ , 
which leads to a significant simplification of the equations. On the other hand the power-law 
structure of the lagrangian make the equations simpler. 

Approximate solutions have been found in three regimes of approximation which describe the 
evolution of the ratio of the energy densities of the electric and magnetic fields during inflation. 
It is assumed that initially the energy density of the electric and magnetic field are of the same 
order. Furthermore, these initial fields are due to quantum fluctuations in the electromagnetic 
field during inflation. Whereas in the radiation dominated era, the energy density in the magnetic 
field decreases as a -4 , the electric field strength rapidly decays in the highly conducting plasma 
(see, e.g., [HES])- Solutions in closed form have been found and the resulting primordial magnetic 
field estimated. It has been shown that depending on the regime of approximation and the value 
of the Pagels-Tomboulis parameter 5 primordial magnetic fields can be generated that are strong 
enough to seed a galactic dynamo. Thus we have provided an example of a theory of nonlinear 
electrodynamics where the nonlinearities act in a way as to amplify sufficiently an initial magnetic 
field. 

The energy-momentum tensor of the electromagnetic field can be cast in the form of an imperfect 
fluid. This has been found explicitly for the particular model of nonlinear electrodynamics under 
consideration here. Moreover, this allows to find the expression for the energy density p of the 
fluid. Requiring that p should be positive provides the bound 5 > \. 

In [29] the possible creation and amplification of magnetic fields was studied in an inflationary 
model coupled to a pseudo Goldstone boson (see also [?)). In this case the lagrangian has the form 
L ~ ^d^Od^e -X + g a 6Y, where 9 is the axion field. This provides an example of a more general 
lagrangian having also an explicit dependence on Y = \F^ U *F tJlU . However, as it turns out the 
resulting primordial magnetic field is not strong enough in order to seed, for example, a galactic 
dynamo. In [30J the model of [29] was generalized to N axions. In this case it was found that at 
least the weaker bound of r > 10 -57 can be satisfied. Here, in this work the creation of primordial 
magnetic fields in a particular model of nonlinear electrodynamics has been studied. It might be 
interesting to generalize this to lagrangians depending on both electromagnetic invariants X and 
Y. 
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6 Appendix A 



In this section it is checked that the approximate exact solution (J376J) is a good approximation 
to the solution of the full differential equation (|3.2p . The solution (|3.6p satisfies equation (|3.5|) . 
Writing the full differential equation (|3.2p as 



yy = 5y 2 + 



m 



(6.1) 



where for the approximate solution y = C2 cosh(z) 1 ~ s with z = m(x + (5 — l)Ci) the additional 
term I is given by 



/ = 



r 28+l 

On 25 + 1 

cosh 1 ~ s (z) 



ai(S- 1) 
4(5 + 1) / a 



(25 - 1) 



m 



x 1 



m 



1 



w n 20 / x 
tanh(z) H — ^ — 



1-5) 2 
-2 



tanh (z) 
2 



+ 



m 
T 



(6.2) 



At X2 when the comoving length scale A leaves the horizon z = by construction. Thus / is 
proportional to ^ 1- At the end of inflation, x = x±, using the bound on —mx\ which in 

general implies —mx\ <C 1, I(x±) is given approximately by 



I(xi) 



20Cf +1 , 2s+i 
— — 2 cosh »-i (zi), 



(6.3) 



where the last factor is much less than 1 since it is assumed that 5 > 1 and, moreover, z\ ~ 
—mx\e N ^ 3> 1. Thus choosing C2 appropriately, |/(xi)| <C 1. 

Finally, it can also be checked using the bounds on — mx\ that the square of the magnetic field 

n 



strength B? is well approximated by equation 



7 Appendix B 

Expressions for j^- and j^- for the solution ()3.30|) . 
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cosh 



-2 



m 



Cl 



(7.1) 



(7.2) 



where /i = — (^) 2 x\. 



19 



References 

[1] P.P. Kronberg, Rep. Prog. Phys. 57 (1994) 325; E. G. Zweibel and C. Heiles, Nature 385 
(1997) 131; L. M. Widrow, Rev. Mod. Phys. 74 (2002) 775; M. Giovannini, Int. J. Mod. Phys. 
D 13 (2004) 391; R. Wielebinski and R. Beck (Eds.) "Cosmic Magnetic Fields", Lect. Notes 
Phys. 664 (Springer, Berlin Heidelberg, 2005); M. Giovannini, "Magnetic fields, strings and 
cosmology," |arXiv: astro- ph /06 1 2378 . 

[2] M.J. Rees, Q.J.R. Astron. Soc. 28 (1987) 197. 

[3] D. Grasso and H. R. Rubinstein, Phys. Rept. 348 (2001) 163; A. D. Dolgov, "Magnetic fields 
in cosmology," arXiv:astro-ph/0306443. 

[4] M. S. Turner and L. M. Widrow, Phys. Rev. D 37 (1988) 2743. 

[5] B. Ratra, Astrophys. J. 391 (1992) LI; D. Lemoine and M. Lemoine, Phys. Rev. D 52 (1995) 
1955; M. Gasperini, M. Giovannini and G. Veneziano, Phys. Rev. Lett. 75 (1995) 3796. 

[6] O. Bertolami and D. F. Mota, Phys. Lett. B 455 (1999) 96; A. Ashoorioon and R. B. Mann, 
Phys. Rev. D 71 (2005) 103509. 

[7] M. Giovannini, Phys. Rev. D 62 (2000) 123505; K. E. Kunze, Phys. Lett. B 623 (2005) 1. 

[8] F. D. Mazzitelli and F. M. Spedalieri, Phys. Rev. D 52 (1995) 6694; C. G. Tsagas and A. Kan- 
dus, Phys. Rev. D 71 (2005) 123506. 

[9] M. Born, Nature 132 (1933) 282; Proc. Roy. Soc. A 143 (1934) 410; M. Born and L. Infeld, 
Nature 132 (1933) 970 and 1004; Proc. Roy. Soc. A 144 (1934) 425. 

[10] W. Heisenberg and H. Euler, Z. Physik 98 (1936) 714. 

[11] J. Schwinger, Phys. Rev. 82 (1951) 664. 

[12] Z. Bialynicka-Birula and I. Bialynicki-Birula, Phys. Rev. D 2 (1970) 2341. 

[13] S. L. Adler, J. N. Bahcall, C. G. Callan and M. N. Rosenbluth, Phys. Rev. Lett. 25, 1061 
(1970). S. L. Adler, Annals Phys. 67 (1971) 599. 

[14] A. K. Harding, M. G. Baring and P. L. Gonthier, Astrophys. J. 476 (1997) 246; M. G. Baring 
and A. K. Harding, Astrophys. J. Lett. 507 (1998) L55; M. G. Baring, Phys. Rev. D 62 (2000) 
016003; A. K. Harding and D. Lai, Rept. Prog. Phys. 69 (2006) 2631. 

[15] E. S. Fradkin and A. A. Tseytlin, Phys. Le tt. B 163 (1985) 123 ; A. A. Tseytlin, "Born-Infeld 
action, supersymmetry and string theory," |arXiv:hep-th/ 9908105 . 

[16] G. W. Gibbons and C. A. R. Herdeiro, Class. Quant. Grav. 18 (2001) 1677; Phys. Rev. D 63 
(2001) 064006. 

[17] R. G. Leigh, Mod. Phys. Lett. A 4 (1989) 2767. 

[18] J. Plebahski, "Lectures on non-linear electrodynamics", (NORDITA, Copenhagen, 1970). 
[19] G. I. Rigopoulos, E. P. S. Shellard and B. J. W. van Tent, Phys. Rev. D 73 (2006) 083521. 



20 



[20] A. A. Starobinsky, In *De Vega, H.j. ( Ed.), Sanchez, N. ( Ed.): Field Theory, Quantum 
Gravity and Strings*, 107-126; A. S. Goncharov and A. D. Linde, Sov. Phys. JETP 65 (1987) 
635; H. E. Kandrup, Phys. Rev. D 39 (1989) 2245. 

[21] H. Arodz, M. Slusarczyk and A. Wereszczynski, Acta Phys. Polon. B 32 (2001) 2155; M. Slusar- 
czyk and A. Wereszczynski, Acta Phys. Polon. B 34 (2003) 2623. 

[22] H. Pagels and E. Tomboulis, Nucl. Phys. B 143 (1978) 485. 

[23] G. F. R. Ellis, "Cargese Lectures in Physics" Vol.6, ed. E. Schatzman (New York: Gordon and 
Breach) pi; C. G. Tsagas and J. D. Barrow, Class. Quant. Grav. 14 (1997) 2539; C. G. Tsagas, 
Class. Quant. Grav. 22 (2005) 393; J. D. Barrow, R. Maartens and C. G. Tsagas, Phys. Rept. 
449 (2007) 131. 

[24] K. Subramanian and J. D. Barrow, Phys. Rev. D 58 (1998) 083502. 
[25] A. Dolgov, Phys. Rev. D 48 (1993) 2499. 

[26] A. C. Davis, M. Lilley and O. Tornkvist, Phys. Rev. D 60 (1999) 021301. 
[27] S. Hannestad, Phys. Rev. D 70 (2004) 043506. 

[28] M. C. Bento, O. Bertolami and N. J. Nunes, Phys. Lett. B 427 (1998) 261; M. C. Bento and 
O. Bertolami, Phys. Lett. B 365 (1996) 59. 

[29] W. D. Garretson, G. B. Field and S. M. Carroll, Phys. Rev. D 46 (1992) 5346. 

[30] M. M. Anber and L. Sorbo, JCAP 0610 (2006) 018. 



21 



